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ABSTRACT 

We measure the angular two-point correlation and angular power spectrum from the NRAO VLA Sky Survey (NVSS) of radio galax¬ 
ies. They are found to be consistent with the best-fit cosmological model from the Planck analysis, and with the redshift distribution 
obtained from the Combined EIS-NVSS Survey Of Radio Sources (CENSORS). Our analysis is based on an optimal estimation of 
the two-point correlation function and makes use of a new mask, that takes into account direction dependent effects of the observa¬ 
tions, sidelobe effects of bright sources and galactic foreground. We also set a flux threshold and take the cosmic radio dipole into 
account. The latter turns out to be an essential step in the analysis. This improved cosmological analysis of the NVSS emphasizes the 
importance of a flux calibration that is robust and stable on large angular scales for future radio continuum surveys. 
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1. Introduction 


Continuum radio galaxy surveys can probe large cosmic vol¬ 
umes. The mean redshift of observed radio galaxies is signifi¬ 
cantly higher than the mean redshift of modern optical surveys. 
Optical surveys allow us to obtain photometric or even spectro¬ 
scopic redshift measurements, which is not the case for contin¬ 
uum radio galaxy surveys. Nevertheless, the information con¬ 
tained in continuum radio galaxy surveys can be useful for prob¬ 
ing cosmological models. An example is the investigation of 
non-Gaussianity via angular two-point correlations. In this work 
we revisit some aspects of the cosmological analysis of contin¬ 
uum radio galaxy surveys. 

Most of our current understanding of cosmology relies 
on observations of the cosmic microwave background (CMB). 
Planck satellite data extended our understanding to the high 
multipole CMB, and allowed us to successfully measure the 
cross-correlation between the large-scale structure and the CMB 
through the integrated Sachs-Wolfe (ISW) and lensing effects. 
In the study of the ISW effect, large-scale radio continuum cata¬ 
logues play an important role. For the Planck analysis the NRAO 
VLA Sky Survey (NVSS) (|Condon et al.|1998|) has been utilized 
dPlanck Collaboration. XIX.|2014|). 


Radio galaxy surveys, such as the NVSS, have large sky 
coverage and extend to redshifts well beyond unity. Their an¬ 
gular two-point correlation function and angular power spec¬ 
trum are useful probes of large-scale cosmology. With new radio 
continuum surveys, containing between 10 million to a few bil¬ 
lion objects, such as from the Low Frequency Array (LOFAR^^ 
the Australian Square Kilometre Array Pathfinder (ASKAP^^ 
the South African MeerKAT radio telescope (MeerKAT^^or the 
Square Kilometre Array (SKAQ more complete and precise ra- 
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dio point source catalogues will be available. These surveys will 
allow us to probe several interesting aspects of cosmology (for a 
recent review see [Jarvis et al.| ( [MT5] )), including accurate inves¬ 


tigations of the largest cosmic structures. 

For cosmological analysis the fundamental observable from 
a survey is the surface density of sources above a given flux den¬ 
sity limit, i.e. the integral of the differential number count per 
solid angle per flux density interval as a function of position on 
the sky. A study that includes all effects of the surface density 
at linear order in cosmological perturbation theory has been pre- 


sented recently (Maartens et al. 

20131 [Chen & Schwarz[[2015[). 

These effects could be modifiec 
tures, e.g. by local voids ([Rubar 
surveys are also a good candidj 

by local effects of large struc- 
t et al.|2014). Radio continuum 
ate to further test several CMB 

anomalies ([Bennett et al.|2011[[Copi et al.|2010[ Naselsky et al. 


2012 


et al. 


[Planck Collaboration XXIII|2014l Copi et al.|2015a Copi 
2015bl [Planck Collaboration. XVI|2015[)~ 


Previously, [Blake & Wall[ ( [2002b[ ) explored the angular two- 
point correlation of the NVSS catalogue at small angular sep¬ 
aration (6 < 10°). In addition, they als o measured the an gular 
power spectrum of the same catalogue ( [Blake et ^[2004[ ), and 
found that their result seems to be incompatible with a high red¬ 
shift population of radio galaxies. [Xia et al] ( [2010[ ) suggested 
that the clustering at 1° < ^ < 8° is dominated by some non- 
Gaussian component. They quote the non-Gaussianity parame¬ 
ter 25 < /ni < 117 at 2c r. The angular power spectrum w as also 
tested more recently in [Marcos-Caballero et al.[ ( [2013[ ), which 
also find an excess of power at large scales, but -43 < /ni <142 
at 2cr. [Hemandez-Monteagudo[ ( [2010[ ) investigated in detail the 
two-point cross-correlation between the NVSS catalogue and the 
CMB temperature anisotropies from the Wilkinson Microwave 
Anisotropy Probe (WMAP) in angular and multipole space. He 
did not find any evidence for the cross-correlation in the range 
I € [2,10] where, according to the simulation 50% of the signal 
is expected. He concluded that the large-scale clustering excess, 
which he found in the NVSS catalogue, is unlikely to be caused 
by contaminants or systematics since it appears to be indepen¬ 
dent of the flux density threshold. 


Article number, page 1 ofpT] 


















































A&A proofs: manuscript no. aa 



0*00 


+175* 


Fig. 1. Surface density of the NVSS source catalogue in galactic coordinates in Mollweide projection, shown at healpix resolution Nside - 32. 
The colour bar shows the surface density in units of number of objects per square degree. Here we include all objects contained in the catalogue. 
The white dots indicate the the positions of the celestial poles. 


In this work, we explore the full range of the angular two- 
point correlation function (0° < ^ < 180°) with a focus on 
large angular scales (> 1°) and compare it with the flat A cold 
dark matter (ACDM) model. The focus of our investigation lies 
on how to prepare a catalogue from a real all-sky radio contin¬ 
uum survey for cosmological analysis, rather than the estima¬ 
tion of cosmological parameters. All previous analyses are based 
on masking of the galax y by a constant latitude cut. Moreover 
the cos mic radio dipole ( [Ellis & Baldwin|p~984t [Blake & Wall 


2002a[ ) has not been accounted for in previous works. To extract 


the auto-correlation signal from the catalogue in a reliable way, 
we propose a more sophisticated NVSS mask and subtract the 
cosmic radio dipole. Based on sidelobe contamination estimates 
and noise properties of the NVSS catalogue, we provide a mask 
that leaves us with a sky coverage of 64.7% at a minimal speciflc 
flux density of 15 mJy. 

Our work is structured as follows. In section [2] we review 
the properties of the NVSS catalogue and explain our masking 
method. In sectionj^ we review the theoretical angular two-point 
correlation function based on the ACDM model. We present the 
angular two-point correlation obtained by means of the Landy- 
Szalay estimator using our mask in section In section [5] we 
discuss our results, followed by the conclusion in section!^ 


2. NVSS catalogue 

The NRAO VLA Sky Survey [Condon et al.[ ( [1998[ ) used the D 
and DnC antenna conflgurations of the Very Large Array (VLA). 
The survey was carried out between 1993 and 1997. Contin¬ 
uum intensity and linear polarization images at 1.4 GHz were 
obtained, covering the whole northern and part of the southern 
sky at declinations S > -40°. The D conflguration of the VLA 
was used to cover the sky from S = -10° to b = 78°, the rest was 
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fllled in by means of the DnC conflguration. The obtained im¬ 
ages have on average a LWHM resolution ^fwhm = 45", which 
is signiflcantly larger than the median angular size of faint ex- 
tragalactic sources (around 10"). This survey design sacriflces 
resolution to achieve high surface brightness, which is needed to 
achieve flux-limit completeness. 

Since the NVSS point-source response is much larger than 
the median angular size of extragalactic sources, most of the 
information on the NVSS total intensity images is represented 
well by elliptical Gaussian fits. The fitted parameters formed 
the NVSS catalogue of sources (not all of which are individual 
sources). The NVSS catalogue contains almost 2 million dis¬ 
crete sources with flux density above 2.5 mJy. Below we use the 
source positions in right ascension a and declination S, the inte¬ 
grated flux density S , and their errors. 

The NV SS catalogue contains several artificial effects ( [Ho[ 
[et al.[[20Q8] ): Lirst, the catalogue shows a conflguration effect, 
which is easily seen in figure as steps in surface density at 
b ~ 78° and S ~ -10°, which are the borders between the D 
and DnC conflgurations on the map. Second, the sidelobes of the 
bright sources can obscure the presence of weak nearby sources. 
They are not completely removed by cleaning. The NVSS local 
dynamic range of the total intensity images is about 1000 to 1. 
Thus, sources closer than about 0.6° to a bright source of flux 
density 5 and fainter than 10“^5 can only be caused by sidelobes 
( [Condon et al.|199^ . 

In addition to these two unphysical effects, there is also a 
foreground of point sources mainly from the Milky Way (see 
figure and smooth components of radiation from the galaxy 
itself. These contaminations can influence the VLA noise tem¬ 
perature, and further change the completeness in certain areas on 
the sky. 
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s > 

5 mJy 

10 mJy 

15 mJy 

25 mJy 


160.1 

59.7 

30.7 

21.5 


Table l.;^f^-values testing for the isotropy in declination of the NVSS 
surface density after masking the region \b\ < 5° and for 13 degrees of 
freedom. 


All these effects need to be treated carefully, otherwise the 
statistical analysis based on the contaminated map will be bi¬ 
ased. Here we choose to apply a lower flux density threshold 
and to mask the catalogue to address these issues. 


2.1. Lower flux density threshold 


Two different VLA configurations (D and DnC) have been used 
to compile the NVSS catalogue. The VLA C configuration is less 
compact and has a resolution of 15", which is better than 45" of 
the D configuration. The DnC configuration is a hybrid configu¬ 
ration in which the antennas on the east and west arms are in D 
configuration, but those on the north arm are in C configuration 
to enhance their view of sources at low elevation. Using the DnC 
configuration changes the synthesized beam and the resolution 
of declination with respect to the D configuration. This shifts the 
brightness sensitivity and completeness between the parts of sky 
observed with different configurations. 

The source catalogue is derived from intensity images, there¬ 
fore it is brightness sensitivity limited. Apparently, the D config¬ 
uration has higher brightness sensitivity than the DnC configu¬ 
ration. Thus a more complete catalogue for the part of the sky 
observed by means of the D configuration is expected. 

The completeness shift of the DnC configuration sky with 
respect to the D configuration sky can be considered as a faint 
source selection based on the position angle, noise, confusion 
and cataloguing procedures. It is not clear whether this selection 
also introduces a tension in source distribution between the two 
parts of the sky. Therefore, it is safe to either use one configu¬ 
ration alone (which would reduce the sky coverage by a signif¬ 
icant amount), or to choose a higher flux density threshold for 
the cosmological analysis (which reduces the source density and 
increases the shot noise). 

[Blake & Wall] ( |2002b| ) argue that this configuration effect is 
only significant at flux densities S < 10 mJy. In figure we plot 
the surface density fluctuation. Act I a = o-ja - 1, as a function 
of declination. The dependence of the surface density fluctuation 
on declination resembles the theoretical rms noise level of the 
NVSS catalogue ( jCondon et al.||1998| ), which is discussed fur¬ 
ther in section!^ It can be seen that the declination dependence 
of the surface density is less than 2.5% and the mean value of 
the DnC configuration is clearly lower than that of the D config¬ 
uration. This configuration effect is more pronounced at lower 
thresholds. 

This finding is also supported by a -analysis testing for 
deviations from isotropy in declination 




V)in 


i=l 


A'bin 

AU_ = yQ. 

(■>!*],r - 


(Ao-,-)^ 


( 1 ) 


where Q/ denotes the solid angle of the ith declination strip. The 
results of this test are shown in table [T] 

In this work, we choose flux density thresholds of 15 mJy 
and 25 mJy for cosmological analysis for the following reasons: 



Fig. 2. Surface density fluctuation, Acr/cr, of the NVSS catalogue as 
a function of declination for several flux density thresholds. The error 
bars assume Poisson distributed source counts and we mask the region 
\b\ < 5°. For clarity, the horizontal positions of the 5 >5 mJy and 
S' > 15 mJy data points are slightly offset. 


The study of the configuration effect suggests that there is a sig¬ 
nificant dependence on declination. As is shown above, a flux 
density threshold of 15 mJy or 25 mJy reduces the value of;^^ 
dramatically compared to the situation with lower thresholds. As 
we do not expect a perfect agreement with isotropy, it is not jus¬ 
tified to raise the threshold even higher. 

Our choice is supported by studies of the cosmic radio dipole 
from NVSS, which is one to two orders of magnitudes larger 
than the higher multipole moments. The previous NVSS dipole 
measurements show that for thresholds below 15 mJy, the dipole 
direction differs significantly from the CMB dipole direction, 
and the reduced of the quadratic dipole estimator increases 
significantly ([Blake & Wall||2002at [Gibelyou & Huterer|[2012t 
[Rubart & Schwarz|2Q13t[Rubart|2015|). 

No cut-off at high flux densities is introduced as there is only 
a small number of sources at high flux densities and thus they 
only play a subdominant role. We also expect them to be less 
affected by calibration systematics. 


2.2. Masking strategy 


Nevertheless, we have to mask regions around the brightest 
sources in order to minimize contamination from their sidelobes. 
In the analysis of [Ho et~ST ( 2008 ) all sources above 2.5 Jy 
have been masked by a disk with a radius of 0.6°. However, not 
all sidelobes appear as spurious entries in the catalogue, a uni¬ 
form masking strategy will erase all the information from bright 
sources, since all sources above the cut threshold are masked. In 
Blake & Wall[ ( [2002b| ) a list of 22 masks around bright galaxies 
was compiled based on a visual inspection of the survey. 

Here, we introduce an automatic bright source selection 
method. We count the number of nearby sources (within 0.6°) 
with 5 > 15 mJy of the brightest sources whose flux densities 
are S > 2.5 Jy. The corresponding histogram is shown in Fig.[^ 

We assume that the number of sources in a fixed solid an¬ 
gle follows a Poisson distribution, an d use the so-called Pois- 
sonness plot ( [Hoaglin & Tukey|2006[ ) to identify the histogram 
bins that significantly deviate from a Poisson distribution, where 
we assume that such a deviation is caused by the sidelobe con¬ 
tamination of bright sources. The idea of the plot is to consider 
a simple variable substitution, which transforms the exponen¬ 
tial Poisson distribution function to a linear function. We And 
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Fig. 3. Cumulative distribution of nearby source counts (disks with 
radius 0.6°). The red line corresponds to nearby source counts around 
bright radio galaxies with S > 2.5 Jy. The blue line corresponds to 
1000 randomly picked positions outside the galactic plane (\b\ > 5°). 
The maximum of randomly picked counts results in 49 nearby sources. 


NVSS position uncertainty map 



Fig. 5. NVSS position uncertainty map, relative to the mean beam 
width ^fwhm “ 45^^. 

15 mJy, the rms position uncertainty is dominated by noise. Ac¬ 
cordingly, we create a position uncertainty map (Fig. 0 by aver¬ 
aging the position uncertainties 



Fig. 4. Poissonness plot of the nearby source counts. Ao is the mean 
nearby count over the survey area. N is the total number of sources 
with S > 2.5 Jy. The solid horizontal line corresponds to a perfectly 
Possion distributed nearby source count. The error bars denote the 99% 
confidence levels. 


that 49 sources fail the Poissonness test at the 99% confidence 
level (Fig. 0. It is worth clarifying that ‘clean’ regions contain¬ 
ing bright sources that, by chance, contain the same amount of 
sources as ’dirty’ regions are also excluded. We also verify that 
for most of the bright sources the cumulative nearby source dis¬ 
tribution is in good agreement with the cumulative distribution 
of randomly picked nearby source counts, which justifies the in¬ 
clusion of many of the bright sources in our analysis. 

We now turn to the issue of noise and confusion. According 


Condon et al. (19981, the rms position uncertainty cr^ 


pos 


IS 


O' b^FWHM 

25 , ’ 


( 2 ) 


where 5p is the peak flux density and (Xb is the rms brightness 
fluctuation (noise and confusion). Our idea is to use CTp to trace 
(Tb- [Condon et al.| ( [l998| ) point out that for flux densities below 
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0-0 = 


I O-^O-q; 

I ^FWHM 



( 3 ) 


for all point sources in a pixel whose flux density is smaller than 
15 mJy. We note that the dominant sources are those with low 
flux density. The map is constructed using the healpix package 
with the pixel size flxed by Nside = 32. 

The resulting position uncertainty map is shown in Fig. 
Which clearly shows that the position uncertainty follows the 
theoretical rms noise level which, outside of the galactic plane, 
increases away from the zenith owing to pickup of ground 
radiation, atmospheric emission, ionospheric effects and uv- 
projection. The galactic plane and nearby nebulas are also easy 
to identify. We also employ a 5% pixel cutoff for the pixels with 
highest 0-0 > 0.132. In addition, we mask galactic radio sources 
by excluding all sources with galactic latitude \h\ < 5°. 

To sum up, at 15 mJy and 25 mJy thresholds we mask all 
pixels in the neighbourhood of 49 selected bright sources. Addi¬ 
tional pixels and their neighborhood with highest mean position 
uncertainty ctq and galactic sources with \h\ < 5° are masked 
as well. In the following we call this the NVSS65 maslQ It is 
shown in Fig.|^ In total, approximately 64.7% of the sky is left. 
The total number of objects after applying the NVSS65 mask at 
15 and 25 mJy is shown in Table 


3. Angular two-point correlation: theoretical 
expectation 


The angular two-point correlation function w{ 6 ) is a powerful 
tool for measuring the projected large-scale structure distribu¬ 
tion of the Universe. It is defined as the joint probability dP of 
finding galaxies in both of t he elements of solid angle and 
6^2 separated by an angle 6 ( Peebles|1980 ), 


dP = d-^dQ.^ 8 D. 2 [\+w{e)]. 


(4) 


^ The NVSS65 mask is public available as a source file on the arXiv 
1507.02160. 
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Fig. 6. NVSS65 mask. The orange region makes up 64.7% of the sky. Yellow pixels are close to the galactic plane \b\ < 5° and are excluded 
to suppress galactic point sources and foregrounds. Blue pixels are excluded owing to large position errors. Red and black pixels contain bright 
sources with significant sidelobe effects; black pixels also overlap with pixels with high position error. 
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Fig. 7. Surface density of the NVSS source catalogue for a flux density threshold of S' >15 mJy and applying the NVSS65 mask, shown in 
galactic coordinates at pixel size Nside-32. The colour bar shows the surface density cr in units of number of objects per square degree. 


At linear order, the full relativistic expression of the surface den¬ 
sity of sources above flux density is ( Chen & Schwarz|2Q15| ), 




Js, 




^phy 


2-r (Of-r 


[1 -i- Sfi + 2if/ + 


. 6r dSr 


+y'^^;* + 2— -r 


dro 


■ 2/Cg 


where So is the observed flux density, a is the scale factor, Tq 
is the inferred comoving distance, a is the mean source spec¬ 
trum index, is the conformal Hubble parameter, and fiphy is 
the physical mean number density. Inside the brackets, the most 
important contribution comes from the number density fluctua- 


(5) 
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tion 6n. The other terms describe volume distortion via the scalar 
metric perturbations j/r and E\ light cone projection involving 
the projected source velocity radial distance fluctuation 

IdrjrQ + ddr/dvo, which lead to a fluctuations of flux density; 
and lensing convergence Kg. We only take the dominant leading 
term, i.e. physical number density fluctuation, into account and 
neglect all other effects. This is justified since the statistical and 
systematic fluctuations from the NVSS catalogue are well above 
the expected size of the subleading linear effects. 

The luminosity and density evolution of radio galaxies is sig¬ 
nificant. Therefore, we turn to redshift space. 


o-(> S t) 


C _ 

' I dz-^nphy[l 


+ 6 „ 


( 6 ) 




dcr 

dz—[\+b6], 

dz 


where b = b{z) denotes the bias of radio galaxies and 6 is the 
matter density contrast in synchronous comoving coordinates. 

We use the Combined EIS-NVSS Survey Of Radio Sources 
(CENSORS) ( [Brookes et al.|2008] ) to model the redshift distribu¬ 
tion of the NVSS catalogue. CENSORS contains all the NVSS 
sources above 7.2 mJy that are within a patch of 6 deg^ in the 
ESO Imaging Survey (EIS). Eollowing [Marcos-Caballero et al.j 
( 2013| ), we choose the gamma function redshift distribution. 





(7) 


with TV denoting a normalization factor that is irrelevant for the 
final result and the best-fit parameters to CENSORS data given 
by zo = 0-53 ^q| 3 and (5 = O.SI^q^^. However, these numbers 
should be treated with care as the redshift distribution is only 
based on 149 galaxies. 

Owing to the broad shape of the luminosity function, the ra¬ 
dio galaxy redshift distribution shows only a weak dependence 
on the flux threshold ( jde Zotti et al.||2Q10l [Blake et ari|2004| ). 
This argument holds true, because for 5 >1 mJy the differential 
number counts are dominated by active galactic nuclei, i.e. the 
flux density threshold is well above the flux density of a typical 
star forming galaxy. 

Eor a statistically isotropic distribution of radio galaxies the 
angular two-point correlation w{6) can be determined from the 
angular power spectrum C/, 

W{0) = T ^(2/ + DQPKcos 9 ), (8) 


where Pi{x) are Legendre polynomials. The angular power spec¬ 
trum Cl of the surface density fluctuation can be calculated from 
the underlying matter power spectrum, 

dk 

—P{k)T{kff{l,k)\ (9) 

with P{k) denoting the primordial power spectrum of density 
fluctuations, T{k) the matter transfer function and 

f(l,k) = J' ^b(z)ji[ck7}(z)]D(z)dz, (10) 

where (dd-/dz) is the differential number density of sources at 
flux density above St, b(z) denotes the bias, j}(x) is a spherical 
Bessel function of first kind, p = ?](z) denotes conformal time, 
and D(z) is the growth factor. 
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Eor the galaxy bias, we use the second-order polynomial 
from the Planck 2015 analysis of the integrated Sachs-Wolfe ef¬ 
fect ( [Planck Collaboration. XXI.|2015[ ), which is an approxima¬ 
tion of the Gaussian bias evolution model of [Xia et ar| ( [2010| ): 


Z7(z) = 0.9 fl+0.54(1+z)2j. 


( 11 ) 


Recently, a detailed analysis of the biasing and evolution of 
NVSS galaxies was presented by [Nusser & Tiwari| ( [2015[ ). Their 
best-fit bias function is comparable to Eg. ([11[) up to z 1.5 and 
grows slower at larger redshifts. We checked that this difference 
is insignificant for the purpose of this study. Eor the add itional 
bias due to non-Gaussianity see[Matarrese et al.[([2000[) and[Dalal[ 

radTO . -— — 


The expected Q for the NVSS catal ogue are obtained us ing a 
modified version of the class package ( [Di Dio et al.|2013| ). The 
parameters for the best-fit cosmological model are taken from 
Planck Collaboration. XVI. | ( [2014[ ). Eor the theoretical predic¬ 
tion of the angular two-point correlation we cut the Legendre 
series at /max = 900. We convinced ourselves that this cut-off 
is large enough to ensure numerical convergence of w(0) at all 
angular scales considered in this work. 


4. Angular two-point correlation of radio galaxies 

4.1. Optimal estimation 


The angular two-point correlation could be estimated from (Pee- 
[bles|1980[ ) 


1 -l- w(^) — 


DD 


Q 


V(V-r l)/2 {Skiy 


( 12 ) 


where DD denotes the count of pairs at separation 6 and TV de¬ 
notes the total number of objects considered in the analysis. Thus 
N(N l)/2 is the total number of possible pairs, kl is the solid 
angle of the survey, and {Ski) is the averaged solid angle of the 
ring 6 to 6SO within kl for a randomly placed ring centre in kl. 

Eor our analysis we use the optimal estimator found by 
[Landy & Szalay[ ( p"993[ ), 


>^LS(^) - 


Vr(TVr -^1)DD 


DR 

■ V(Vr + 1)= + 1, 

RR 


(13) 


where RR means the averaged pair count over a number of large 
random simulations, DR is the averaged data-random cross pair 
count for a number of large random simulations, and denotes 
the number of sources in the random catalogues. It is necessary 
to clarify that the simulated random sources R have to be identi¬ 
cal for RR and DR to minimize the statistical uncertainty. Under 
the assumption that higher order correlations among galaxies can 
be ignored, the Landy-Szalay estimator has a “Poisson error” 


dwLs(^) - 


1 1-1- wls(^) 

^[DD 1 + wq 


(14) 


where wq is the mean of the two-point correlation function over 
the sampling geometry. 


Wo 


= r Gp{6)w{6)dkl 

Jq 


(15) 


where Gp{6) is a dimensionless geometric form factor which is 
equal to the fraction of unique c ell pairs separated by distance 
6 ± do 12 ( [Landy & Szalay|199^ . 
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N 

d 

a 

S 

S > 10 mJy 

5 > 15 mJy 

S > 25 mJy 

436,733 

314,594 

200,092 

1.32 X 10 -^ 
1.44 X 10-2 
1.83 X 10-2 

142.70^ 

153.44° 

157.12° 

30.47° 

-5.53° 

-15.10° 

expected 


0.46 X 10-2 

o^ 

00 

0 

-7° 


Table 2. NVSS dipole for various flux density thresholds, measured by 
means of healpix at resolution Nside - 32 after applying the NVSS65 
mask. For comparison we quote the expected kinetic dipole for NVSS 
radio sources, based on the observed CMB dipole. 


In |Xia et ar] ( |2010| ) a pixel based estimator was used {Nside = 
64). We compared the Landy-Szalay estimator and the pixel 
based estimator at Nside = 64 for the same flux threshold and 
mask. On angular separations above 5 degrees, we found good 
agreement between both estimators. However, for smaller angu¬ 
lar scales it seems that the pixel estimator is redistributing cor¬ 
relation from the smallest scales to intermediate scales. As for 
Nside = 64 the distance between two neighbor pixel centers is 
as large as 0.9 deg, it is not surprising that a pixellation error 
shows up at the few degrees scale. 


4.2. Dipole correction 

The radio dipole signal is believed to result from our peculiar 
motion ( [Ellis & Baldwin [19 84 ) with respect to the cosmic rest 
frame of radio galaxies, because that the mean redshift of ra¬ 
dio galaxies is above one. If this rest frame is the same as the 
CMB rest frame, then the dipole measured in the radio catalogue 
should agree with the CMB dipole measured by WMAP and 
Planck dJarosik etal. 120111 IPlanck Collaboratio m XXVII. 20141). 
Previous studiej^ ^ Blake & Wall|2002at|Singalf20lH|Gibelyou| 


& Huterer j2012[ Rubart & SchWrz||26l3t [Tiwari et al.||2015t 

Singal|2014| ) measured the radio dipole for the NVSS catalogue. 
It is actually signiflcantly larger than expected, by a factor of 
two to four depending on the details of the analysis. The origin 
of this dipole excess is currently unknown. One possibility might 
be that it is a combination of local large-scale structure ( jRubartj 
et al.||201^ and a kinetic component due to Doppler shift and 
aberration. 

A local structure dipole, and the kinetic dipole as well, vio¬ 
late the assumption of statistical isotropy that is implicit in the 
way we estimate w(6). We suggest that the dipole signal needs 
to be taken into account prior to the further correlation or power 
spectrum analysis. In our analysis the largest contribution to the 
angular two-point correlation of radio galaxies also comes from 
the dipole moment of the galaxy distribution, as can be seen in 
Fig.[^ A signiflcant effect can also be observed for the extracted 
multipole moments at low I (not shown here). 

Utilizing healpix, we And the radio dipole of the NVSS cat¬ 
alogue after masking with NVSS65 (see table [^. The estimated 
dipole at our chosen flux density thresholds an d sky coverage 
agrees with the estimates from the literature (see jRubartj ( j2015j ) 
for a recent summary). 

The standard dipole subtraction approach for pixelized maps 
is relatively straightforward. First, we estimate the dipole ampli¬ 
tude and direction through a linear dipole estimator on the pix¬ 
elized map. Then we subtract the measured dipole contribution 
at each pixel. However, we do not use the pixel map to measure 
the correlation function, but rather extract it using the measured 


^ Note that the definition of the radio dipole amplitude in Blake & Wall 
( 2002aj ) differs from the rest of the literature by a factor of 2. 


xIO"" 



Fig. 8. NVSS angular two-point correlations for a 5° galactic latitude 
mask with and without dipole correction. 


positions of all radio sources outside the mask. To achieve this, 
we have to include the effect of a dipole into the Landy-Szalay 
estimator. 

Thus we simulate random catalogues with the measured 
dipole, denoted by below, and employ the following estima¬ 
tor. 




NriNr +\)DD 

N{N+\) m 


- N(N, + 1)^ + 
RR 




(16) 


The dipole simulations Rd are achieved by the following proce¬ 
dure. First, we assign a uniform random number t to each sim¬ 
ulated object and then modify this number based on the angular 
separation if/ between the object and dipole direction. 


t = Random[0,1) + 2 ^ 


(17) 


Then we add the objects with t > 0.5 to the catalogue and 
drop the others. Each random catalogue contains N^ = 10^ ob¬ 
jects and we average over ten such catalogues. 

The dipole modified estimator ( p^ makes use of the full po¬ 
sition information of the sources and by simulating several large 
random catalogues, we minimize the uncertainty in w^^(6). The 
computational load is the disadvantage of this procedure. 

We can now compare the results with and without dipole 
subtraction. We either corrected for the measured radio dipole 
(see table 1^ or for the CMB predicted radio dipole. The dipole 
contribution in the angular two-point correlation function can be 
seen from Fig. We find that the dipole has a significant ef¬ 
fect and actually dominates the two-point correlation function at 
large angular scales above ~ 10°. In the figure we account for 
the measured dipole. Considering just the CMB predicted dipole 
reduces the large-angle correlation, but leaves us with a residual 
dipole that could be due to a local structure and that is hard to 
predict without a much more detailed study. We therefore de¬ 
cided to correct for the measured NVSS radio dipole and also 
suppress the structure dipole in the theoretical prediction. 


4.3. Results 

We adopted the dipole subtracted Landy-Szalay estimator to 
measure the angular two-point correlations of the NVSS cat- 
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\b\ < 5“ 

NVSS65 


N 


N 


S > 15 mJy 

377,739 

165.96 

322,557 

94.08 

S > 25 mJy 

240,872 

222.49 

205,103 

99.80 


Table 3.;^^-test for w{6) for 49 data points [excluding first bin {0 < 6 < 
3.6°)]. We neglect the correlation between the data points. 


alogue with thresholds of 15 mJy and 25 mJy for two differ¬ 
ent masks (NVSS65 and a constant latitude cut of the galactic 
plane). 

The results agree with our expectation. After eliminating the 
contribution from high rms noise and confusion pixels by means 
of the NVSS65 mask, the two-point correlation turns out to be 
less scattered, (see Fig.|^. For the simpler constant latitude cut, 
the first data point at ~ 4° is above the plot range of the figure. 
The NVSS65 mask efficiently reduces the amount of correlation 
at scales of a few degrees and brings the measurement in agree¬ 
ment with the theoretical expectation of the best-fit cosmological 
model, and also observe a better agreement of the 15 and 25 mJy 
thresholded data sets after the NVSS65 mask has been applied. 
These findings are confirmed by the;^^-values shown in Tablej^ 
We infer that the new NVSS65 mask efficiently pushes the data 
points towards the theoretical prediction. 

Fig. also shows that the measurement using the NVSS65 
mask agrees with the theoretical prediction until 6 ~ 90°. Above 
that angular scale the data appear to be more noisy. The quite 
large correlations at the largest angular scales close to 180° are 
surprising. 

For an angular scale between 1° and 20°, the 65% sky cov¬ 
erage guarantees that the w(6) estimator is averaged over a 
large number of independent sky patches. Artificial fluctuations 
caused by the survey or galactic foreground are suppressed in 
the average, and sidelobes and multicomponent source effects 
are expected to contaminate smaller angular scales (up to 0.6°). 
In this range, w(6) is expected to be consistent with the theoreti¬ 
cal prediction. 

The result of an analysis at higher angular resolution for 
6 < 20° is shown in Fig. [T^ In order to suppress the shot noise 
contribution we now focus on the 5 > 15 mJy data set. We And 
that the NVSS65 mask improves the agreement with theoretical 
predictions considerably. At this angular scale, the most impor¬ 
tant effect of the NVSS65 mask is to remove spurious correlation 
at ^ < 2 deg. We also observed that the NVSS65 mask increases 
the importance of the dipole subtraction at a separations of a few 
degrees. As jXia et al.| ( |201()| ) apply a flux threshold of 10 mJy 
in their analysis, we also show our corresponding analysis for 
comparison, however as already discussed above, we think that 
this flux threshold is not a conservative choice. A discussion of 
the cosmological consequences is given below. 

A complementary analysis to the angular two-point correla¬ 
tion function is to study the angular power spectrum. One way of 
measuring C/ is to do a Legendre transformation of the angular 
two-point correlation. The result of that transformation is shown 
in Fig. m Many of the C/ turn out to have negative values, which 
shows that this measurement is quite noisy. Nevertheless, in the 
mean the C/ seem to agree well with the theoretical expectations, 
with exception of a few multipoles with even /, most prominently 
the / = 10 mode. When restricting the Legendre transformation 
to the range ^ < 100 deg, assuming w(0) = 0 at 6 > 100 deg, we 
And that the most extreme values, e.g. at ^ = 10 are suppressed. 
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Fig. 9. Angular two-point correlation function w(6) from 5° galactic 
latitude cut (top panel) and the NVSS65 mask (bottom panel). In both 
cases we include a dipole correction. 



Fig. 10. Angular two-point correlation function at 1° < ^ < 20° from 
the NVSS catalogue with dipole correction for two different masks and 
two different flux thresholds. 


5. Discussion 

5.1. Cosmological implications 

We now turn to the cosmological implications and a discus¬ 
sion of our findings in comparison with previous studies. In 
the following we focus on the results obtained by means of 
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Fig. 11. NVSS angular power spectrum C/ for S' > 15 mJy, dipole 
corrected and NVSS65 masked. C/ is evaluated via a Legendre trans¬ 
formation from the angular two-point correlation function (red). The 
extreme fluctuations are reduced, if the Legendre transformation is lim¬ 
ited to ^ < 100 deg (green). The solid line and the band around it show 
the theoretical prediction and its cosmic variance. 



Fig. 12. Angular two-point correlation function for different redshift 
distributions. 


the NVSS65 mask, including the dipole correction as discussed 
above and a lower flux threshold of 15 mJy. 

A question of interest is the determination of the cosmologi¬ 
cal parameters of our minimal six-parameter cosmological stan¬ 
dard model. For that purpose, the NVSS data alone cannot com¬ 
pete with high fidelity data from the CMB or from optical galaxy 
redshift surveys. Therefore it is not interesting to determine any 
of the six cosmological base parameters from a fit to NVSS. 

It is important to note that the redshift distribution function 
also affects the NVSS angular two-point correlation at this angu¬ 
lar scale. The Figure [T^ shows that the parameter zo changes its 
slope and and angular scale. Our angular two-point correlation 
measuremen t seems to prefer low zo and yg , which agrees with 
the results of |Marcos-Caballero et al.| ( [2Q13| ). 

There is a clear degeneracy between cosmological parame¬ 
ters such as Hq, the redshift distribution parameters (zo and JS) 
and the galaxy bias, which limits our ability to use the most re¬ 
cent radio continuum surveys for cosmological parameter esti¬ 
mation. However, it is encouraging that the completely indepen¬ 
dent measurements of Hq (Planck), zo and {3 (CENSORS) pro¬ 
vide a picture that seems to be consistent with w{6) from NVSS. 
In order to extract the full science potential of upcoming radio 
surveys, a larger sample of redshifts must be measured. 

To go beyond the set of the six base parameters, the small an¬ 
gle two-point correlation is also an interesting probe of primor¬ 
dial non-Gaussianity. Previous results ( |Xia et aH20I0| ) claimed 


evidence of a primordial non-Gaussianity at the 3cr level. The 
angular correlation function is supposed to vanish at around 4° 
if the fluctuations are Gaussian, but their measurement of w{6) 
showed a constant shift from zero at 1° < ^ < 10°, which they 
attributed to the effect of a primordial non-Gaussianity. How¬ 
ever, adapting our procedure of masking, dipole correction and 
optimal estimation of w{6), we do not see this shift. Our result 
is consistent with primordial Gaussianity, as is shown in Fig.p^ 
Turning that into a precise limits on f^\ is beyond the scope of 
this work. 


A primordial non-Gaussianity would also lead to an increase 
of angular power C/ at low multipoles ( |Xia et aL||20I0t |Des 


[Jacques & Seljak||20I0[ ). The angular power spectrum obtained 
via a Legendre transformation is shown in Fig. m One could be 
tempted to conclude t hat we do find an in crea se of power for low 
multipoles, similar to [Blake et ar| ( |2004| ) and [Marcos-Caball^ 


[et al.[ ( [20T3] ). However, it should be noted that our method to es¬ 
timate the angular power spectrum is not optimal and appears 
to be noisy. On top, most of the excess power is in parity even 
modes, primordial (local) non-Gaussianity would not care about 
parity. 


We tried to understand what causes for the removal of the 
apparent non-Gaussianity in the correlation function, but it turns 
out that it is impossible to attribute it to a single effect. As al¬ 
ready discussed in detail above, the essential steps are, using 
the Landy-Szalay estimator (to avoid pixellation errors), using 
a larger mask, and removing the dipole contribution. 


5.2. Residual systematics 

The theoretical prediction based on the ACDM model suggested 
Cl ^ 5 X 10“^. Our measurements closely agree with that pre¬ 
diction, with some exceptions. The quadrupole cannot be de¬ 
tected at any significant level and the power at / = 10 is an 
order of magnitude larger than expected. However, at larger /, 
up to / ~ 60, the Cl are consistent with the theoretical prediction, 
which is in agreement with previous analyses of the ISW effect 
Hernandez-Monteagudo[ ( [2QT0| ). Higher multipole moments are 
noisy and statistically consistent with zero. 

We finally discuss the angular scales at ^ > 20° and look in 
more detail at the corresponding multipole moments up to / = 
10. For a simple galactic isolatitude cut, we find significantly 
more power at low multipole moments, with I = 4 being the 
dominant mode (not discussed here). The essential way to get 
rid of this extra power is to take direction dependent systematic 
effects into account in the NVSS catalogue. Our NVSS65 mask 
allows us to reduce the / = 4 mode and to recover an overall 
fiat angular power spectrum. However, even with the NVSS65 
mask, some of the Qs are one order of magnitude larger than 
the prediction at / < 10. Perhaps because of remaining surface 
density fluctuations at different declinations. 

We find a clear anti-correlation between the surface density 
and the theoretical rms noise of the NVSS catalogue. Fig, 

The declination dependence of the theoretical rms noise fluctu¬ 
ations at -40° < S < -10° is due to changes in the snapshot 
integration time. When the VLA points to the horizon, the ef¬ 
fect of the projection of the uv plane combined with ground and 
ionospheric noise decreases the effective signal-to-noise ratio of 
the survey. As can be seen in Fig. this effect is not limited 
to the faintest sources, but is also there for brighter sources at 
5 > 60 mJy. Thus, only the most extreme influence of the di¬ 
rection dependent systematics can be cured by means of a lower 
flux threshold. 
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Fig. 13. Top panel: Surface density fluctuation (Acr/^) of NVSS 
sources at A >60 mJy and A > 15 mJy after applying the NVSS65 
mask. Bottom panel: The theoretical rms noise level of the NVSS 
catalogue from [Condon et al.| (|1998|). The VLA geodetic latitude is 
34°04'44". 


To further probe whether these direction dependent effects 
affect the low-/ multipoles, we investigated the angular pseudo¬ 
power spectra, which are obtained through a decomposition into 
spherical harmonics on the pixelized galaxy number count. For 
this we used healpix. In figure we show estimates of the 
pseudo-C/ for two different normalizations. In the first case, we 
normalize the pixel count with the mean pixel count of the full 
survey n, i.e. 


Xi 


Hi - n 


n 


( 18 ) 


The second case is motivated by |Vielva et al.| ( [2006| ). We divide 
the map along the declination, based on durations of the snap¬ 
shots, and subtract the mean pixel count hs^ at declination A/, 



Fig. 14. The healpix pseudo-C/ of the NVSS catalogue with NVSS65 
mask for A > 15 mJy and radio dipole subtracted. The zig-zac line 
shows the standard normalization. The crosses show the declination nor¬ 
malization result. The smooth solid line shows the theoretical prediction 
and the band around it includes shot noise and cosmic variance. 



e (degree) 


Fig. 15. Pixel based estimate of the angular correlation function of the 
NVSS catalogue for A > 10 mJy with \b\ < 5 mask (red line). The ef¬ 
fects of dipole removal (green) and declination normalisation (brown) 
are shown for comparison. The purple line shows the result using a dec¬ 
lination normalization and subsequent dipole subtraction. 




■ns, 


This procedure was also used in the Planck analysis of the ISW 
effect (|Planck Collaboration. XIX.|2014[ [Planck Collaboration.] 
[XXI.|20l5l ); 

This second procedure, erases the source density fluctuations 
along the declination direction. As can be seen in Fig. 

/ = 3,4,5 multipoles are significantly reduced by the dec 
mean subtraction, which strongly implies that the I 


14 the 


Imation 
4 mode 

fluctuations are caused by the direction dependent noise of the 
NVSS catalogue. On the other hand, the large / = 10 mode is not 
affected by that normalization at all and it is not clear why it ex- 


(19) pixels) as a baseline in Fig. 15 This baseline seems to have some 


excess power. We found that both the strip normalization and the 
dipole subtraction give a substantial reduction and almost iden¬ 
tical results. We thus think that it would be better to just sub¬ 
tract the dipole and to use a more conservative mask or carefully 
modeling this systematic effect, instead of setting all fluctuations 
between declination strips to zero. 

In addition, we note a curious observation. Looking carefully 
at Fig. m we observe a parity asymmetry in the C/: the power 
of even multipoles / = 4, 8,10,12,16,22,24,32... is higher than 
that of their odd neighbours. This even parity preference can also 
be seen because w(180°) > 0. A similar effect was observed 
in the CMB angular power spectrum, where the parity asym- 


nation normalization for a cosmological analysis, such as[Vielva 

metry is the opposite Kim & Naselsky] ([2010|). There the odd 

et al. (2006[), 

Fernandez-Cobos et al.[([2014[) or [Planck Collab- 

/(/ -f l)C//27rs are larger than the even ones. 

oration. XIX. 

(2014[), certainly suppresses fluctuations on large 

Our results further suggest that w(6) from full-sky radio con¬ 
tinuum surveys can be used to constrain cosmological parame- 

scales, but the procedure cannot distinguish direction dependent 


effects from real fluctuations and it is very hard to assign an error 
estimate to it. 

To analyse this effect in more detail, we took the A > 10 
mJy, \b\ < 5 degree map before dipole subtraction and used the 
same pixel estimator as in [Xia et al.] ( [2010[ ) (including only full 


ters at an epoch that is barely accessible to other probes. Planned 
and upcoming surveys with instruments like LOFAR, ASKAP, 
MeerKAT and Anally SKA will allow us to reduce the shot noise 
and to increase angular resolution and sensitivity, while cover¬ 
ing all sky, extending the studies to several frequency bands (our 
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discussion here is limited to 1.4 GHz), and improving the con¬ 
trol of systematic effects. For more details see [Raccanelli et alT 
( |2012| ) and [Jarvis et al.| ( |20T5] ) and references therein. 


6. Conclusion 


In this paper we revisited the angular two-point correlation func¬ 
tion w(6) and angular power spectrum C/ from the NVSS cata¬ 
logue of radio galaxies. The focus of our work was to investigate 
systematic effects in the NVSS catalogue. In order to minimize 
the contribution of these effects in the cosmological analysis, we 
provide a new NVSS mask with 64.7% percent of sky, called 
NVSS65. We also find that it is essential to account for the radio 
dipole and to use an optimal estimator. We found that our mask 
significantly improves the;^^ value of the angular two-point cor¬ 
relation function on all angular scales. For angular scales be¬ 
tween 1 and 20 degrees, w(6) agrees with the flat ACDM model 
without introducing primordial non-Gaussianity. 

Thus we have shown that to fully explore the cosmological 
potential of continuum radio surveys, one has to understand and 
investigate the systematic effects related to flux calibration, espe¬ 
cially direction dependent effects of the calibration. In addition, 
the effect of the cosmic radio dipole affects the reconstruction of 
higher multipole moments and the attempts to measure or con¬ 
strain primordial non-Gaussianity. 

To obtain an improved upper limit on f^\ or to constrain other 
cosmological parameters at a redshift of about unity is beyond 
the scope of this work, as it would need an extensive study of 
the uncertainties coming from our understanding of the density, 
luminosity and bias evolution of radio galaxies. However, our 
analysis shows that the radio sky is in remarkably good agree¬ 
ment with the standard model of cosmology according to Planck. 
It will also be interesting to improve the ISW analysis of the 
Planck-NVSS cross-correlation by means of the new NVSS65 
mask, and to include a radio dipole correction, as well as a higher 
flux threshold. 


Acknowledgements. We thank Anne-Sophie Balleier, Daniel Boriero, Bin Hu, 
Dragan Huterer, Hans-Rainer Klockner, Alvise Raccanelli, and Matthias Rubart 
for valuable comments and discussions. We acknowledge financial support from 
Deutsche Forschungsgemeinschaft (DFG) via the Research Training Group 
“Models of Gravity” (RTG 1620). We are grateful for the possibility of per¬ 
forming the numerical computations on the Nordrhein-Westfalen state comput- 
ing cluster at RWTH Aachen. We acknowledge the use of the NVSS catalogue 
[Condon et al. (|199^, provided by the NR AO. This work made extensive us e 
of the HEALPi x Gorski et'n!]j2005) and the class [Lesgourgues & Tram[{2011b 
|Di Dio et al.H2013^ packages. 


Ellis, G. F. R. & Baldwin, J. E. 1984, MNRAS, 206, 377 
Fernandez-Cobos, R., Vielva, R, Pietrobon, D., et al. 2014, MNRAS, 441, 2392 
Gibelyou, C. & Huterer, D. 2012, MNRAS, 427, 1994 
Gorsld, K., Hivon, E., Banday, A., et al. 2005, ApJ, 622, 759 
Hernandez-Monteagudo, C. 2010, A&A, 520, AlOl 

Ho, S., Hirata, C., Padmanabhan, N., Seljak, U., & Bahcall, N. 2008, 
Phys. Rev. D, 78, 043519 

Hoaglin, D. C. & Tukey, J. W. 2006, Exploring Data Tables, Trends, and Shapes 
(John Wiley & Sons, Inc) 

Jarosik, N., Bennett, C., Dunkley, J., et al. 2011, ApJS, 192, 14 
Jarvis, M., Bacon, D., Blake, C., et al. 2015, PoS, AASKA14, 018 
Kim, J. & Naselsky, P. 2010, ApJ, 714, L265 
Landy, S. D. & Szalay, A. S. 1993, ApJ, 412, 64 

Lesgourgues, J. & Tram, T. 2011, J. Cosmology Astropart. Phys., 1109, 032 
Maartens, R., Zhao, G.-B., Bacon, D., Koyama, K., & Raccanelli, A. 2013, J. 
Cosmology Astropart. Phys., 1302, 044 

Marcos-Caballero, A., Vielva, R, Martinez-Gonzalez, E., et al. 2013, ArXiv e- 
prints 1312.0530 

Matarrese, S., Verde, L., & Jimenez, R. 2000, ApJ, 541, 10 
Naselsky, R, Zhao, W, Kim, J., & Chen, S. 2012, ApJ, 749, 31 
Nusser, A. & Tiwari, P. 2015, ArXiv e-prints 

Peebles, P. 1980, The Large-Scale Structure of the Universe (Princeton Univ. 
Press) 

Planck Collaboration. XIX. 2014, A&A, 571, A19 

Planck Collaboration. XVI. 2014, A&A, 571, A16 

Planck Collaboration. XVI. 2015, ArXiv e-prints 1506.07135 

Planck Collaboration. XXL 2015, ArXiv e-prints 1502.01595 

Planck Collaboration XXIII. 2014, A&A, 571, A23 

Planck Collaboration. XXVII. 2014, A&A, 571, All 

Raccanelli, A., Zhao, G.-B., Bacon, D. J., et al. 2012, MNRAS, 424, 801 

Rubart, M. 2015, PhD thesis. University Bielefeld 

Rubart, M., Bacon, D., & Schwarz, D. J. 2014, A&A, 565, Alll 

Rubart, M. & Schwarz, D. J. 2013, A&A, 555, A117 

Singal, A. K. 2011, ApJ, 742, L23 

Singal, A. K. 2014, A&A, 568, A63 

Tiwari, R, Kothari, R., Naskar, A., Nadkarni-Ghosh, S., & Jain, P. 2015, As- 
troparticle Physics, 61,1 

Vielva, P, Martinez-Gonzalez, E., & Tucci, M. 2006, MNRAS, 365, 891 
Xia, J.-Q., Viel, M., Baccigalupi, C., et al. 2010, ApJ, 717, L17 


References 

Bennett, C. L., Hill, R. S., Hinshaw, G., et al. 2011, ApJS, 192, 17 
Blake, C., Ferreira, P. G., & Borrill, J. 2004, MNRAS, 351, 923 
Blake, C. & Wall, J. 2002a, Nature, 416, 150 
Blake, C. & Wall, J. 2002b, MNRAS, 329, L37 

Brookes, M., Best, P, Peacock, J., Rottgering, H., & Dunlop, J. 2008, MNRAS, 
385,1297 

Chen, S. & Schwarz, D. J. 2015, Phys. Rev. D, 91, 043507 
Condon, J. J., Cotton, W. D., Greisen, E. W., et al. 1998, ApJ, 115, 1693 
Copi, C. J., Huterer, D., Schwarz, D. J., & Starkman, G. D. 2010, Adv. Astron., 
2010, 78 

Copi, C. J., Huterer, D., Schwarz, D. J., & Starkman, G. D. 2015a, MNRAS, 449, 
3458 

Copi, C. J., Huterer, D., Schwarz, D. J., & Starkman, G. D. 2015b, MNRAS, 
451,2978 

Dalai, N., Dore, O., Huterer, D., & Shirokov, A. 2008, Phys. Rev. D, 77, 123514 
de Zotti, G., Massardi, M., Negrello, M., & Wall, J. 2010, A&A Rev., 18, 1 
Desjacques, V. & Seljak, U. 2010, Classical and Quantum Gravity, 27, 124011 
Di Dio, E., Montanari, F., Lesgourgues, J., & Durrer, R. 2013, J. Cosmology 
Astropart. Phys., 1311, 044 


Article number, page 11 ofpT] 
















